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ABSTRACT 



We present estimates of the N-point galaxy, area-averaged, angular correla- 
tion functions uiNi^) for N = 2,. ..,7 for galaxies from the fifth data release of 
the Sloan Digital Sky Survey. Our parent sample is selected from galaxies with 
^ ■ 18 < r < 21, and is the largest ever used to study higher-order correlations. We 

subdivide this parent sample into two volume limited samples using photometric 
redshifts, and these two samples are further subdivided by magnitude, redshift, 
^ and color (producing early- and late-type galaxy samples) to determine the de- 

^ • pendence of ujn{G) on luminosity, redshift, and galaxy-type. We measure u}n{G) 

■ using oversampling techniques and use them to calculate the projected, sat. Using 
CNl ■ models derived from theoretical power-spectra and perturbation theory, we mea- 
t:J" ' sure the bias parameters hi and C2, finding that the large differences in both bias 
^ ■ parameters {hi and C2) between early- and late- type galaxies are robust against 
O ■ changes in redshift, luminosity, and as, and that both terms are consistently 

smaller for late-type galaxies. By directly comparing their higher-order correla- 
^ ■ tion measurements, we find large differences in the clustering of late-type galaxies 

■ at redshifts lower than 0.3 and those at redshifts higher than 0.3, both at large 
scales (c2 is larger by ~ 0.5 sX z > 0.3) and small scales (large amphtudes are 
measured at small scales only for z > 0.3, suggesting much more merger driven 
star formation at z > 0.3). Finally, our measurements of C2 suggest both that 
a 8 < 0.8 and C2 is negative. 

Subject headings: Cosmology: observational; Galaxies: clustering, environment, 
star formation 
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Introduction 



An important subject in cosmology and galaxy formation is galaxy bias. Galaxies trace 
the underlying dark matter that shapes the Universe, but there is no guarantee that they 
do so faithfully. The galaxy bias defines the relationship between the clustering of dark 
matter and the clustering of galaxies. Ideally, the relationship between the distribution of 
galaxies and dark matter is simple — the ratio of the over-density of galaxies to the over- 
density of dark matter (i.e., the bias, hereafter b, which is measured relative to clustering 
dark matter) is a constant factor independent of smoothing scale and the over-density of 
dark matter. This ideal case is known as scale invariant linear bias. If instead 6 is a function 
of the overdensity, it is non-linear. Significantly non-linear bias would suggest that galaxy 
formation is dependent on environment, as it would suggest that the efficiency of galaxy 
formation is dependent on halo mass (i.e., the dark matter over-density). 

Assuming a non-linear bias, b can formally be expanded into factors b^ via a Taylor 
expansion, which reduces to linear bias whe n b^ = b and b^ = for N > 1. Several previoius 



studies have found that 62 is nearly zero. |Pan fc Szapudil (120051 ) found, by studying the 
monopole contribution to the 2dFGRS three-point c orrelation function, that, for galaxies 
with -21 < bj < -20, 62 = -0.06l[J:|][5?. Likewise, Iflikage et all J2005h found, using the 
Sloan Digital Sky Survey (SDSS) spectroscopic galaxy catalog to calculate the bispectrum, 
that C2 (cat = bj\f/bi) is consistent with zero to within 10% for erg = 0.9, ho wever, they 



measure C2 to be larger for smaller cxg- On the other hand, in an earlier study we (iRoss et al. 



20061 . hereafter denoted Ross06) found C2 = —0.30 ±0.10 via determination of the Nth-order 
angular galaxy correlations for 11 million photometrically selected galaxies from the third 
data release (DR3) of the Sloan Digital Sky Survey (SDSS). iGaztafiaga et al.l (120051 ) likewise 
found C2 = — 0.36^Q Qg for bj < 19.45 galaxies in the 2dFGRS by calculating the redshift-space 
three-point correlation functions, a method that is independent of ag, the rms fluctuation in 
the matter density averaged in a sphere of radius 8 Mpc, . 



Determining cig is necessary, because, as yet, there is no clear c onsensus. When me asured 
from only the third year WMAP (WMAP3) data, erg = 0.744+°;°^° JSpergel et al.lbooeh . Fur- 
ther constraining WMAP3 data via inflationary models, cxg = O.TOOlg oes. When combined 
with the SDSS power-spectrum data, however, this increases to 0.772j^ nnj?s, and wh en com- 



bined with the 2dFGRS power-spectrum data it decreases to 0.737]'^o;o45 (ISpergel et al.ll2006[). 
Furth ermore, the SDSS power spectrum by itself determines Ug 



0.842+°:°^^ (iTegmark et al.l 



2OO4J ). M ethods that employ clus t ers of galaxies typically measure smaller values of erg. For 
example, IVoevodkin fc VikhlininI (120041 ) found o"g = 0.72 ± 0.04 by using a cluster baryon 
mass function. 



Independent measurements of bi and ag are not possible using correlation measurements. 
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as as serves as the normalization to the model power spectra one must use to determine bi. 
The value of 62? however, is not degenerate with ag, as it will change the shape of correlation 
measurements. A novel technique to provide an independent measure of ag is to capitalize on 
the independence of 62- By determining both 62 and C2, Ug can be constrained. Determining 
the value of C2 alone places loose constraints on the value of cxg. 

Our most dramatic findings in Ross06 were that significant non-linear bias was required 
to explain the observed clustering of early- and late-type galaxies. The results determined 
that early-ty pe galaxies' Co was la rger than that of late- type galaxies by an absolute difference 
of 1.0±0.13. ICroton et al.l (120061 ) have since confirmed this measurement — that red galaxies 
have a higher C2, via their determination of higher-order correlations in the 2dFGRS, though 
they found a more modest difference of 0.36±0.17. In another study, iNishimichi et al.l (120061 ) 
found a difference of closer to 0.5 when studying the bispectrum of galaxies selected from 
the SDSS data release 4 (DR4) spectroscopic survey. It is not surprising that differences 
are observed between these three studies, as each uses a different color cut to separate the 
red and blue samples . The existence of non-linear bias between galaxy types was first found 
by iGaztanagal (1l992l ) , when it was discovered that the higher-order clustering of galaxies 
selected by optical surveys (CfA and SSRS) differs from the clustering of galaxies selected 
by the infrared survey IRAS (galaxy selection in IRAS was biased towards late- type galaxies) . 
This was confirmed by lFry fc Gaztanagal (Il993l ) and interpreted as IRAS (and thus late-type) 
galaxies having a negative C2 relative to optically selected galaxies. 

Clearly, galaxy type affects the non-linear nature of the bias. The linear bias is also 
strongly dependent on galax y type, as red, early-type, galaxies have consistently been shown 



to cluster more strongly (e.g., Willrner et al.lll998 : Norberg et al.ll2002l : iMadgwick et al.ll2003 



Zehavi et al.ll2005l : IRoss et al.lboodTlCroton et al.H2006h t han blue, late-type, galaxies. Lumi 



nosity has been shown to scale proportionally to 61 (e.g.. IMadgwick et al.ll2003l : iZehavi et al. 



20051 ) . which logically follows if one expects more luminous objects to generally be more 
massive. It is also likely that the properties of galaxies, and thus b, evolve over cosmic 
time. Therefore, a true characterization of galaxy bias, both linear and non-linear, requires 
quantifying any dependence on galaxy type, luminosity, and redshift. 

As a result, we employ the wealth of information available from SDSS DR5 to follow-up 
and vastly improve our DR3 measurements on the nature of galaxy bias. Apart from Ross06, 
four other s tudies h ave used SPSS phot ometr ic data to calculate h igher-order correlations: 
Gaztanagal (l2002al ) , iGaztaiiagal (l2002bl ) , and ISzapudi et al.l (120021 ) measured higher-order 
correlations for galaxies in the SDSS Early Data Release (EDR) and found that measure- 
ments using S DSS were consistent with previous results and were free of any systematics. 
Blaizot et al.l (120061 ) measured the DRl higher-order correlations and found good agreement 
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with simulation. The RossOG measurements employed the largest data set ever used to cal- 
culate higher-order correlations with over 11 million galaxies from DR3. The SDSS DR5 
offers an opportunity to significantly improve the DR3 results, as it offers ~ 70% more data 
and accurate photometric redshift catalogs complete with rest-frame absolute-magnitudes. 
This wealth of data enables a quantification of the dependence of the Nth-order correlations 
of galaxies on luminosity, type, and redshift. 

In this paper, we, therefore, calculate and analyze the area-averaged angular A^-point 
correlation functions using SDSS DR5 galaxies, up to seventh-order. Our methodology is 
explained in §3. This present work offers a significant improvement over Ross06, due to 
the fact that we use ~ 70% more data, extend the range of measurement by an order of 
magnitude, make measurements using volume limited samples, and employ improved the- 
oretical modeling. Our main sample of galaxies is used for comparison with the Ross06 
measurements, and is split into the same five subsamples as in RossOG, a process which is 
described in §2. The measurements made using DR5 are presented in §4, where we illustrate 
the superiority of DR5 over DR3. We use the SDSS DR5 PhotoZ table to create two vol- 
ume limited samples, allowing us to investigate the evolution and luminosity dependence of 
Nth-order correlations in §5. Model, area- averaged correlation functions are calculated by 
integrating over theoretical power spectra and an appropriate redshift distribution, allowing 
for the calculation of first- and second-order bias parameters. In §6, we present our measure- 
ments of the higher-order (non- linear) bias and determine their dependence on galaxy- type, 
luminosity, and redshift. 

We adopt the cosmology {Qm, ^a, h, T) = (0.28, 0.72, 0.7, 0.15), where T is the shape 
parameter, based on recent supernovae, large-scale stru cture, and CMB measurements (e.g., 
Riess et al.l |2004J : ISpergel et al.l l2006l : I Cole et al.l 120051). This value of the sh ape parameter 
arises naturally (see, e.g.. Equ ations 30 and 31 of lEisenstein fc Hul (119981 )) for a baryon 



fraction of = 0.185 (e.g., ICole et al.ll2005h . 



2. Data 



The data analyzed herein were taken fr om the SDSS DR5 (lAbazajian et al.ll2005l ). This 
surv ey obtains wide-field C CD photometry (ICunn et al.lll998l ) in five passbands {u, g,r,i, z; 



e.g., iFukugita et al.l Il996l ). The entire DR5 represents 8000 square degrees of observing 
area. We selected galaxies with positions lying in the Northern, contiguous portion of the 
SDSS from the DR5 PhotoPrimary database, and further constrained the sample (using the 
Schlegel. Finkbeiner fc David Il998l dust maps) to have reddening-corrected magnitudes in 
the range 18 < r < 21. Further, significant masking was required to account for bright 
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stars and areas of high reddening and poor seeing (see §3.5p . This produced a set of over 
18 milhon galaxies (18,532,911) at a median redshift of about 0.31. This is by far the most 
galaxies used to conduct this type of measurement, and represents over a 70% increase 
in the number of objects used in our previous DR3 measurements. We split this sample 
of galaxies into the same five sub-samples as in Ross06 (three samples constrained to the 
magnitude ranges 18 < r < 19, 19 < r < 20, and 20 < r < 21, and two color-selected 
samples defined by m — r > 2.2 for early- type and u — r < 2.2 for late- type galaxies). These 
samples are used primarily for comparison with the DR3 measurements, however, due to our 
improved measurement techniques, they also sample substantially smaller and larger scales 
than Ross06. 



2.1. Creation of Volume Limited Samples 

Our primary analys is is divided between t wo volume limited samples created following 



the methods outlined in iBudavari et al.l (120031 ). Galaxies with 18 < r < 21 are taken from 



the DR5 PhotoZ table and matched to galaxies in the DR5 PhotoPrimary table. Using the 
rest-frame absolute r-band magnitudes, M^, given for each galaxy in the PhotoZ table, we 
display the Mr - z plane in Figure [TJ As can be seen, there is a definite locus that defines 
the limiting absolute magnitude for a given redshift. To create a volume limited catalog to 
a given redshift, zl, we simply select galaxies with z < zl that are also intrinsically brighter 
than the limiting magnitude defined by the locus of points displayed in Figure [H We focus 
our analysis on two volume limited samples that have roughly similar numbers of galaxies. 
One sample we take is limited to a redshift of 0.4, requiring that < —20.5. This sample, 
hereafter denoted as zA, contains nearly three and half million galaxies (3,380,553) after 
masking (see §3.51) . Our other volume limited sample, hereafter denoted z3, is limited to a 
redshift of 0.3, with Mr < —19.5, and contains nearly four million objects (3,980,652) after 
masking. 



3. Methodology 

3.1. Angular Correlation functions 

We estimate N-point area-averaged angular correlation functions, U]y{9), using a counts- 
in-cells technique (e.g., Ross06). This basically involves calculating the statistical moments 
of the over-densities contained in equal-area cells. The over-density for cell i is defined as 

ft = ^ (1) 

n 
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where n is the average number of galaxies in a cell and is the number of galaxies in cell 
i. The remaining details and equations required to determine u!n{(^) are found in Ross06. 



In a hierarchical model (e.g., iGroth fc Peebles! 119771 : ISzapudi et al.l Il992l : iGaztanaga 



1994 ). higher-order correlations can be expressed in terms of the two-point correlation func- 
tion, and the volume-averaged correlations are given by 

UR) = SnUR)^-' (2) 

where Sn is the hierarchical amplitude. In a similar manner, we can define the analogous 
relationship for the area-averaged angular correlations 



Sn = -r^ 



(3) 



The hierarchical amplitudes of the higher-order moments encode much of the pertinent 
information on the distribution of the data. These amplitudes, therefore, embody the central 
analysis of this paper. 



3.2. Redshift Distributions 



In order to compare angular measurements to theoretical models, it is necessary to 
determine dn/dz, which we accomplish using the galaxies' photometric redshifts (see §2.1). 
We construct dn/dz by using each published redshift and its error (rejecting any with error 
greater than twenty percent) to create a probability density function (PDF). The PDFs for 
each redshift are combined to produce the expected number of objects (n) in a redshift bin of 
width 0.001. This distribution is then normalized and interpolated over in order to estimate 

For SDSS galaxies with 18 < r < 21, the resulting normalized n{z) is plotted in Figure 
El The distribution of n{z) is smooth and roughly Gaussian. For all subsamples (volume 
limited and otherwise, see §2), we likewise use the corresponding DR5 photometric redshifts 
to obtain estimates of their individual redshift selection functions. 



3.3. Pixelization 

The pixelization schemes we employ are nearly identical to those used in Ross06. Basi- 
cally, we reimplemented the SDSSpix pixelization scheme originally developed by Tegmark, 
Xu, and Scrantoiil], as described in detail in RossOG. Increased computing resources have 



^http: / /lahmu.phyast.pitt.edu/^scranton/SDSSPix/ 
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allowed efficient usage of SDSSpix at smaller scales, allowing us to make measurements for 
9 > 0.02°. The striped method (see Ross06) is again employed, but only for scales between 
0.005 and 0.04 degrees (we note that we are able to probe smaller scales than in RossOG). 



For larger angular scales an oversampling technique (e.g.. ISzapudi et al.ll2002l ) is applied 
such that every angular scale uses the same number of data cells. Using the base resolution 
pixels produced by SDSSpix, we can make cells equivalent to any angular scale nOi, at every 
single base pixel (where n is any integer greater than or equal to 2 and 9h is the angular scale 
at the base resolution). Thus, at a large scale the cells are highly overlapping, allowing more 
information to be extracted, which allows for more precise calculations at large scales. We 
therefore perform calculations for 6 < 20°. We have verified that the results using the striped 
method are consistent with using this oversampling implementation at scales between 0.02 
and 0.1 degrees. From here on, we will refer to the method used for small angular scales as 
the striped method and to the method for large angular scales as the oversampling method. 
Oversampling increases the covariance, but this is not a problem since we perform a full 
covariance analysis for all parameters we attempt to measure. 



3.4. Errors and Covariance 



W e compute errors and covariance matrices using a jackkni fe method (e . g..lScranton et al 
2002), with inve rse- variance weigh ting for both errors (e.g., Myers et al. 20051 . 20061 ) and 



covariance (e.g., iMyers et al.l 120071 ). nearly identical to the method described in detail in 
Ross06. The jackknife method works by creating many subsamples of the entire data set, 
each with a small part of the total area removed. For the striped method, we utilize the 
natural geometry of the SDSS. Each of the 29 different stripes in the DR5 forms a natural 
subset of the overall data. The covariance matrix is calculated using each of the possible 
subsamples of DR5 that is made up of 28 stripes. For the angular scales that are calculated 
by the striped method, we find that these 29 subsamples are sufficient to create a stable co- 
variance matrix. For the larger angular scales probed by the oversampling method, we find 
that 20 jack-knife subsamplings are sufficient to create a stable covariance matrix. These 
20 subsamples are created by simply eliminating a contiguous grouping of l/20th of the 
unmasked pixels in 20 separate areas. To properly constrain fit parameters, we minimize the 
using our covariance matrixes via 



X 



(4) 



where C is the covariance matrix and i and j refer to the i*'^ and j*'^ jackknife subsample. 
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3.5. Masks 

We will generally refer to useful observational information (such as seeing and Galactic 
extinction values) across each pixel in our schema (see §3.31) as forming a mask of that 
information. The DR5 area required significant masking which we performed in the same 
manner as in Ross06. Pixels at the base resolution are discarded if they intersect the standard 
SDSS imaging mask, have a mean reddening Ar > 0.2, have a mean seeing greater than l."5, 
or intersect the Ross06 mask for galaxy MlOl. Cells at scales above the base resolution have 
their over-densities corrected for the fractional area of the pixel in the same manner as in 
RossOG, thus 

= (5) 
n 

where Aj is the fractional area of cell i. As before, we did not find any systematic variation 
in correlation measurements as a function of stripe. 



4. Area-averaged Correlation Functions and Hierarchical Amplitudes: 

Complete Sample 

Figure O shows the area- averaged correlation functions for N < 7 determined for galaxies 
with 18 < r < 21. Errors for each point were determined by the jackknife method (see 
RossOG). Compared to the RossOG results, we now sample at more angular scales, and 
the maximum and minimum angles are extended, resulting in an extra order of magnitude 
in angular coverage. For each A^, the correlation function has a shape roughly consistent 
with a power-law to about 0.3°. At larger scales, there exist obvious features that are not 
adequately represented by a power law. Looking at 0)2, a power-law representation appears 
as if it may be appropriate to 1°. Assuming a form A6^~'^ , the best fit for N = 2 over 
the angular range 0.02 < 6 < 1° is A = (7.2 ± 0.1) x 10" ^ and 7 = 1-770 ± 0.005. These 
resul t s are consistent with previous mea surements (see, e.g.,|Gaztanag 3l994l : IConnoUy et al.l 



2OO2I : iFrith. Outram. fc Shanks I l2005bl . RossOG). The value is 433. G, however, meaning 
that a power-law form is inappropriate over these scales. In order to produce fits that are 
not rejected, fits must be performed over ranges that are less than half a dex. We caution, 
therefore that our power-law measurement is useful only for comparison purposes. 

The hierarchical amplitudes si\i{6) for N < 7 measured in magnitude ranges 18 < r < 19, 
19 < r < 20, 20 < r < 21, and 18 < r < 21 are shown in Figure HI Compared to the results 
of RossOG, the measurements made at scales greater than one degree are far more significant, 
but overall the measurements are quite consistent. We also probe smaller scales, showing 
that at scales less than 0.02° the sn measured for 18 < r < 19 increase as the angular scales 
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decreases. 



As in RossOG, we employ the simple color criteria determined by lStrateva et al.l (120011 ) to 
separate early-type {u — r> 2.2) and late-type {u — r< 2.2) galaxies from SDSS photometric 
data. Using this color cut and the magnitude restriction 18 < r < 21, we separate our 
sample into early- and late-type galaxies and repeat our N-point measurements for these two 
samples. We find that DR5 contains about 25% more late-type galaxies (10,284,575) than 
early-type galaxies (8,248,336), which is nearly the same proportions as in Ross06. 

Figures [5] and [6] show the results of the and sn measurements for early- and late- 
type galaxies. The early-type galaxie s clearly show stronger clustering at scales 6* < 5°, in 



agreement with previous res ults (e.g.. IWillmer et al.lll998l : IZehavi et al.ll2002l : iNorberg et al. 



2OO2I : iMadgwick et al.l 120031 ) . but at large scales {9 > 5°), the amplitudes are nearly equal. 
This suggests that their respective bias might be scale dependent (see §6.11) . In Ross06 we 
were not able to probe such large scales, otherwise our overall results remain consistent. The 
correlation functions display roughly power-law behavior for all N, but there is significant 
structure that becomes more pronounced as grows. The of the late-type galaxies show 
smaller amplitudes for 6 > 0.1^, suggesting that higher-order bias terms are significant. At 
the largest scales {6 > 2°), S3 and S4 appear nearly identical for early- and late-types, yet in 
this regime, the errors begin to dominate the measurements. 

As found in Ross06, the late-type correlation measurements show extremely interesting 
behavior. At angular scales between 0.01° and 0.1*^ the late-type galaxies exhibit slopes that 
increase as the scale decreases. Probing smaller scales, this relationship appears to turn over, 
as evidenced by the sn measured at 6 < 0.02°. For A^ = 6 and 7, there is a dramatic loss of 
signal at angular scales greater than 0.2°. A similar loss in signal is shared by the correlation 
measurements made by using all galaxies, but the measurements made by using early-type 
galaxies do not show this dramatic loss of signal. This suggests that there is an intrinsic 
property of late-type galaxy clustering that is so strong it dominates the measurement for 
A^ > 6. 

The DR5 PhotoZ table contains estimates of spectral type, given by a number between 
and 1, where is most red and 1 is most blue. It has been shown (e.g., iBudavari et al. 



20031 ) that splitting the sample of galaxies into early- and late-type galaxies at a type value 
of 0.3 (hereafter deno ted the photo sample) is roughly equivalent to splitting them using the 
Strateva et al.l (I2OOII ) color criteria (hereafter denoted the color sample). To test this, we 
split the entire sample into early- and late-type galaxies by photo and compare the resulting 
correlation functions to those measured for galaxies split by color. The results for S3 are 
displayed in Figure [71 While the two cuts produce similar results, there are some significant 
differences. It appears that the characteristic rise at small scales for late-type galaxies is 
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actually more prevalent in the galaxies split by the photo method. This suggests that type 
confusion when splitting the galaxies by the color method (which would be expected to be 
greatest at higher redshift) may dampen the overall signal of late-type galaxies selected this 
way. To minimize any type confusion, therefore, when we split our galaxy sample by redshift, 
their types will be determined via the photo method, as these types should be less sensitive 
to redshift than the color method (see §5.31 for further justification). 



5. Higher-Order Correlations in Volume Limited Samples 

5.1. z<0.3 

The volume limited sample z3 contains roughly four million galaxies with M,. < —19.5. 
We split this sample into fifteen different subsamples (though each is not mutually exclusive). 
The sample is split by magnitude into two groups: —20.5 < Mr < —19.5 and —21.5 < Mr < 
—20.5, and by redshift into two groups: < ^ < 0.2 and 0.2 < z < 0.3. Each group itself 
contains three groups: early-type, late-type, and all-types, where the split is done by using 
the photo method. 

Figure [8] displays the sn measurements made using all-, late-, and early-type galaxies in 
the z3 volume limited sample. Immediately, one notices that the late-type galaxies do not 
resemble the late-type galaxy measurements made on the full sample, especially at smaller 
scales. The late-type galaxies also display larger error bars than the early-type galaxies, 
despite the fact that there are nearly the same number of late-type (1,984,021) galaxies and 
early-type galaxies (1,996,631) in this sample. The early-type galaxies show significantly 
higher amplitudes than the late-type galaxies, and the difference is largest at ~0.1°, just as 
in the full sample of galaxies. 

For each subset of galaxies (all, late, and early) the amplitudes for the full sample are 
slightly larger than those of z3. This suggests that the bias is slightly larger in the z3 
sample than in the full sample, as sat oc This is not surprising, as the full sample 

conta ins less luminous galaxies than z3, and bias is known to increase with luminosity 
(e.g.. 



Madgwick et al.ll2003l : IZehavi et al.ll2005l ). The bias measurements for each sample are 



calculated in the following section. 

Figure [9] displays the S3 and S4 measurements made by separating the z3 galaxies by 
luminosity into two groups, —21.5 < Mr < —20.5 and —20.5 < Mr < —19.5. These 
measurements for all-, late-, and early-type galaxies are shown in the bottom, middle and 
top panels, respectively. For each galaxy type, the shapes of sn for brighter galaxies are quite 
similar to the fainter galaxies. The largest difference is that the amplitudes are higher for the 
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fainter galaxies, which is consistent with the first-order bias increasing with luminosity. Based 
solely on Figured it does not appear that non-linear clustering (i.e., C2, C3) is dependent on 
luminosity, as the differences in the amplitudes appear consistent with a linear bias model 
(see §6, Tabled]). 

In order to test for evolution, z3 is split by redshift into a.0 < z < 0.2 and 0.2 < z < 0.3 
sample. Figure [10] shows this split for all-, late-, and early-type galaxies in the bottom, 
middle and top panels, respectively. After accounting for the fact that that the galaxies with 
0.2 < z < 0.3 probe scales that are about 1.6 times larger (based on median redshifts of 0.26 
and 0.15) than the < z < 0.2 galaxies, the shapes of the sn are quite similar to a physical 
scale of about 5 /i"^ Mpc (~0.06° for < z < 0.2, ~0.04 for z < 0.3). Despite the visual 
differences, there does not appear to be significant evolution in the bias, as the differences 
in the measurements can be explained by the differences in the physical scales. 

One concern with the volume limited samples we create is the possible effects of pho- 
tometric redshift errors on the creation of our samples. To quantify any potential bias, we 
created ten separate samples where instead of taking the stated photometric redshift of each 
galaxy, we sampled its probability-density-function (PDF) given by its one a error, thereby 
assigning a new redshift for each galaxy. This allowed us to created ten separate samples 
with the same magnitude and redshift limits as z3. Figure [TT] displays the average measured 
S3 of these samples, with errors (in bold black) calculated by finding the standard deviation 
of the ten measurements. Underneath these points, the jackknife errors of the z3 sample are 
plotted in light red. It is clear that the normal jackknife errors dominate the error budget, 
especially at larger scales where the errors are most important for fitting bias values. There- 
fore, we conclude that photometric redshift errors should not make any significant difference 
in our results. 



5.2. z < 0.4 

The volume limited sample zA contains just under three and a half million galaxies 
with Mr < —20.5 and z < 0.4. Figure [12] shows the sn measurements for all-, late-, and 
early-type galaxies in z4, split using the photo method. Their shapes are quite similar to the 
measurements made on the full sample. The z4 sample, however, shows a strong rise at small- 
scales in Sat for late-types, unlike the z3 sample. The error bars on the late-type galaxies 
are again larger that those of the early-type galaxies, but in this case it can be explained by 
the fact that there are significantly fewer late-type galaxies (1,325,488 late-types; 2,055,065 
early-types). 
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Splitting z4 into two redshift bins, z < 0.3 (1,302,750 galaxies) and 0.3 < z < 0.4 
(2,077,803 galaxies) produces extremely interesting results. Figure [13] shows the measured 
in these redshift bins for all-, late-, and early-type galaxies in the bottom, middle, and top 
panels, respectively. For each galaxy-type, the are much closer to constant at the lower 
redshift. For the late- type galaxies, the errors are much larger at small scales for the low 
redshift bin. This is due to the fact that there are nearly 75% more high-redshift late-type 
(843,527) galaxies than low redshift (481,961). Low redshift, late-type galaxies have a weak 
signal at small scales, while the high redshift galaxies have a relatively strong signal, meaning 
that, in the full sample, the high redshift signal dominates. This explains the small-scale 
features seen in the full z4 sample of late-type galaxies. 

The late-type galaxies can be further split by their type. Separating the galaxies with 
type values greater than 0.3 at a type value of 0.65 produces a sample of galaxies that 
corresponds roughly to late-type spirals (0.3 < t < 0.65; denoted LI hereafte r) and one that 



corre sponds roughly to irregular galaxies (0.65 < t; denoted L2 hereafter; iBudavari et al. 



20031 ). Figure [T3] shows the measured S3 and S4 for LI (closed symbols) and L2 (open 
symbols) for 0.3 < z < 0.4. The rise in the amplitudes at small scales is much stronger for 
the L2 galaxies, which we interpret as evidence that the small scale rise correlates with star 
formation in dense environments, which we discuss further in §7.3[ 



5.3. Comparison of Full Sample to Volume Limited 

By comparing the s^v measurements made on the full sample to those made on the 
volume limited samples, features in the full sample of galaxies are isolated in the red- 
shift/luminosity plane. The most obvious conclusion is that the measurements of late-type 
galaxies are dominated by galaxies with z > 0.3, especially on angular scales smaller than 
half a degree, as the sn signal is much stronger at z > 0.3 for late-types. For early-types, 
the shapes of measured in all subsamples are consistent both with each other and with 
the measurement of early- type galaxies drawn from the full sample. 

For many of the subsamples, there appears to be a feature at ~2'^. This is especially 
true for all galaxy types in the main sample of galaxies and for early-type galaxies with 
0.3 < z < 0.4 and z < 0.2. It is quite prevalent for all galaxy types with z < 0.4 and 
is strangely not present for galaxies with 0.2 < z < 0.3. A close inspection of the feature 
reveals that it actually occurs at slightly different angular scales in the different redshift 
shells. Figure [15] shows S3 measured for early-type galaxies in the z < 0.2 redshift bin from 
z3 (red circles) and the 0.3 < z < 0.4 redshift bin of zA (black triangles). These are the 
ideal samples to compare, as they represent the largest difference in physical scale. The left 
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panel shows the two measurements plotted on an angular scale, while the right panel plots 
S3 against the equivalent physical scale (based on their median redshifts of 0.154 and 0.355). 
When the two measurements are plotted against their equivalent scale, the match is much 
better. This suggests that the feature is physical in nature and is characterized by a minima 
at 10 h~^Mpc. While this scale is interesting as it marks the transition from the non-linear 
to linear regime, the error bars preclude us from drawing any profound conclusions. 

Figure [7] suggested that separating early- and late-type galaxies via the color method 
may break down at higher redshifts. Figure [12] displays the same information as Figure 
[T3l with galaxy-type determined via the color method instead of the photo method. The 
measurements for z < 0.3 are nearly identical to those measured splitting by the photo 
method, but for 0.3 < z < 0.4 the results are quite different, for both the early- and late- 
type galaxies. This suggests that for z > 0.3, the color method cannot adequately distinguish 
between early- and late-type galaxies, and justifies our preference for separating galaxy type 
based on the photo method. 



Bias Measurements 



There is no guarantee that galaxies cluster with the same amplitude as the dark matter 
they trace. Indeed, it has often been shown that g alaxies display different clustering am- 



plitudes when separated by luminosity or type (e.g.. iFrv fc Gaztanaga 



19981 : iNorberg et aPbooi iMadgwick etaDl2003l : IZehavi et al" 



2005 



1993 



Ross et al. 



Willmer et al. 



2006 



Croton et al, 



20061 ). The simplest model for this difference is the linear bias model. In this model, the 
overdensity of galaxies is a linear function of the overdensity of dark matter: 



6r, = b6 



DM 



(6) 



where g denotes galaxy, DM denotes dark matter, and b is the linear bias factor. In this 
approximation, one finds the simple relationship that 



DM- 



(7) 



It ha s been shown, howev er (e.g.. lFry fc Gaztanagalll993l : iGaztafiaga et al.ll2005l : IRoss et al. 



20061 : ICroton et al.l l2006l ). that linear bias may not be a good approximation. One can 



represent the relationship between overdensities more generally, such that the measured 
overdensity is some function of t he dark matter overdens ity. As such, the relationship can 
be expanded into a Taylor series (IFry fc Gaztanagal 119931 ): 



00 , 

EbN_^N 



N=0 



DM 
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where ^at is the A^th order bias term. Thus it follows that if non-linear bias is important, a 
measurement of the linear bias will increase as the overdensities increase (i.e., as the scale 
gets smaller) and the non-linear terms (i.e., any for > 1) grow in importance. When 
5 <^ 1, 61 can be determined via Equation [TJ as higher order terms will be negligible. 



6.1. First-Order Bias Measurements 

In order to find bi, one must first determine the scales at which higher-order terms 
make a negligible contribution to U2- This can be done by calculating U2, incorporating a 
second-order bias term, and comparing it to the 0)2 calculated normally. To second order, 
the dark matter over-density can be expressed (via trivial manipulation of Equation [8]) as 



6dm = b^' \^-h ± ^Jbl + 2b26gj (9) 

Thus, when calculating LO2, if Equation [9] is used to correct the measured overdensities in each 
cell, a measurement is returned that assumes a first and second order bias. If the first-order 
bias is set to 1 and a reasonable second-order bias is applied, the resulting measurement will 
begin to deviate significantly from the standard measure ment when the second order effects 



become important. Based on previous measurements (iGaztanaga et al.l l2005l : iRoss et al. 



20061 ) . we select -0.3 as a reasonable value to use. Figure [T71 shows the ratio of UJ2 calculated 
with this value to UI2 calculated in the standard way, using the sample of galaxies from z3 
with —20.5 < Mr < —19.5. The ratio grows significantly greater than 1 for 6 < 0.66°, 
corresponding to a physical scale of about 8/i~^Mpc. This is right where one would expect 
non-linear effects to become important, as it marks the generally accepted transition to the 
weakly linear regime. As a result, we calculate the first order bias for scales where the 
corresponding physical scale (f) is greater than 8/i^^Mpc. 

Calculating the bias requires knowing u}2,dm- Of course thi s is no t directly measurable. 



so we must resort to using a theoretical model. ISmith et al.l (120031 ) have derived fitting 
formulae based on A^-body simulations, which produce matter power-spectra when given 
specific input cosmological parameters. By using a modified version of Limber's equation 



(jPeebleslll980l ). one can use the appropriate redshift distribution to invert the P{k) to obtain 
0O2{e): 

= vr / f ^1 F{x)dz [ P{k)WUDek)dk (10) 



\dz J \dx, 

where W2D = 2^^^^ is the top-hat two-dimensional window function, D is the survey depth 
(determined by the median redshift), P{k) is the matter power spectra, k is the spectral 
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index, and Ji is the first-order Bessel function of tlie first kind. In a fiat universe, F{x) = 1, 
and dz/dx = H{z)/c = -f^oA/^m(l + zY + ^^a/c, simplifying Equation 10 to: 



J (^)V^m(l + ^)3 + f^Arf2 J P{k)Wlj,{Dek)dk (11) 



Tliis integral was numerically determined with our assumed cosmology and fixing the value 
of (jg for the matter P{k) to 0.8. Following this procedure, model uj2 were produced for each 
subsample studied, using the appropriate redshift distribution. Finally, the first order bias 
was calculated by using the covariance matrix at scales f > 8/i~^Mpc. We note that this 
approach is much more accurate than the methods employed in Ross06, which inverted the 
measurements to real-space, required 0)2 be a power-law, calculated only the relative bias, 
and were insensitive to any changes in the bias as a function of scale. 

Tables [U and [2] display our calculated bi values for all of the measurements made using 
the volume-limited samples 2:8 and zA. The most important information can be summarized 
as follows: (1) for z3, there is no significant evolution in hi as a function of redshift for 
any galaxy type, but there is significant evolution seen in z4: and it is most dramatic for 
the late-type galaxies; (2) hi grows larger with luminosity, independent of galaxy type; and 
(3) hi is consistently larger for early-type galaxies, though the ratio of hi^^ariy to &i,/ate is 
larger for z < 0.3. The determined luminosit y and type dependence s of hi are consistent 



with the general results of previous findings (IFrv fc Gaztanagal Il993l: IWillmer et al 



Norberg et al.ll2002l : iMadgwick et al.ll2003l : IZehavi et al 



2005 



Ross et al. 



2006 



1998 



Croton et al. 



20061 ). All of the determined values are dependent on the true value of erg. Since all of the 
models use = 0.8, the true values of the hi are 0.86i,m/c"85 where hi^m are the measured 
values reported in Tables [1] and [2l 

Most of the fits to 61 using the z3 sample minimize such that it is approximately 
one or smaller per degree of freedom (DOF), meaning most fits favor a scale invariant hi 
over the fitted range. The notable exception is the measurement for late-type galaxies with 
0.2 < z < 0.3, with a per DOF value of nearly two. The fits to hi for zA are in general 
worse than those of z3. All but two have fits with x"^ per DOF greater than one. It is thus 
unlikely that hi is scale invariant at higher redshifts, a fact which must be taken into account 
when calculating higher-order bias terms. 



6.2. Second-Order Bias Measurements 



In order to measure the second order bias, we use the equation (IFry fc Gaztanagalll993l ): 

S3,T = K}ri.S3,DM + 3C2,t) (12) 
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where C2 



b2 



This equation was derived for real-space hierarchical amplitudes, but given a 



theoretical ss^dm, C2 can be determined in the same way from g.^ g. Using PT, an e xpression 
valid in the weakly non-linear regime can be derived for Cos^dm ( lBernardeadll995h : 



oo3,DM = 6 J (f )' [n^l + + ^a] dz X 



f (jkP{k) [wio{Dek]f dk) + JkP{k) {w^ij{Dek))dk j k'^Dep{k)W2D{Dek)w'^ij{Dek)dk 



(13) 

From Equation 3, it is clear that s^^dm can be obtained by dividing this uj^^dm by ^jJ2dm- 
Equation 13 uses PT, which means that a linear power spectrum must be used when calcu- 
lating S3,£)Af, and it is only valid at scales with r > 8 /i^^Mpc. Using the model S3, da/ and 
hi calculated at each scale (allowing for any scale dependence in 61), one can calculate C2 at 
each scale and thus construct its covariance matrix and determine the best-fit average 
values (see §3.41) of C2 for each subsample. Once again, we note that this method is superior 



to the one employed by Ross06, as it has been shown (e.g. iGaztanaga fc Bernardeaulll998 
Bernardeau et al.ll2002l ) that this method is a better match to simulations than assuming a 
hierarchy and inverting sa? to real-space. 

The values of C2 for the z3 and zA samples are presented in the right hand columns 
of Tables [1] and [21 The important aspects are as follows: (1) as first seen in Ross06, the 
C2 values for late-type galaxies are significantly lower than those of early-type galaxies; (2) 
early-type galaxies have a value of C2 that is independent of both luminosity and redshift; 
and (3) C2^iate varies slightly with luminosity and evolves significantly between redshifts of 
0.3 and 0.4, which we explore in §7.31 The per DOF for all but three measurements are 
less than 1.24. The three exceptions are for measurements drawn from samples of early-type 
galaxies, and the poor fits are largely due to the strength in the feature at 10 /i~^Mpc. To 
the limit on our data, one would expect most measurements can be fit by a constant C2 , as 
the errors on the S3 measurements over the fit ranges are typically large. 



7. Discussion 
7.1. C2 vs. cTg 

The relationship between C2 and ag is more complicated than for the first-order bias. 
Given that our measured C2,m are calculated using erg = 0.8 for the matter P(/c), it can be 
shown that 
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where C2 is the true c^. Since S3 is not constant, the difference between C2 and C2,m is not 
constant either. Thus, it is possible to perform a two parameter fit on the values of C2 and 
CTg, thereby placing loose constraints on cxg. Figure (TS] shows the one a allowed regions of 
parameter space for C2 and ag for all- (black), early- (red), and late- type (blue) galaxies in 
the z3 volume limited sample. Clearly, there is a large spread in the allowed values of C2 even 
if one constrains 0.7 < cxg < 0.9 (changing ug from 0.7 to 0.9 decreases C2 by about 0.6 for all 
galaxies type). It is clear that a significant difference in the C2 of early- and late- type galaxies 
is robust against changes in cxg, as increasing cig has only a slight effect, (the difference is 1.0 
for cTg = 0.7 and it is 1.1 for erg = 0.9). Thus, irrespective of ag, there is clear evidence for 
non-linear bias differences between early- and late-type galaxies. 

Focusing on the z3 sample of all galaxy types, if we set C2 = we constrain the value of ag 
to be in the range 0.64+0.04/-0.03. This is inconsistent at the one a level with the WMAP3 
best-fit parameters (erg = 0.74;!^Q;Qg) but consistent with the one a lower bound on WMAP3 
as constrained by inflationary models (e.g, ug = 0.700lo;Qg5). On the other hand, it is highly 
at odds with the typical effect of combining WMAP3 and large-scale structure constraints, 
which tend to higher values of ag . We would contend that these higher values of cxg may be 
inconsistent with hnear theory, as our data suggest C2, for all galaxies, devi ates by more than 



two a from zero for ag > 0.71. If we instead set C2 = —0.36, as found by iGaztanaga et al. 



(l2005l ). ag = 0.77 ± 0.04 — nearly identical to the WMAP3-SDSS combined measurement. 
Thus, considering that it is both unlikely that C2 is greatly negative and that it is unlikely 
that ag is less than 0.68, we contend that the best interpretation of our results is that C2 is 
at least slightly negative and that ag < 0.8. 

With the methods employed in this paper, it is not possible to break the degeneracy 
between the bias parameters and ag. It may be possible, however, using measurements 
similar to those presented herein. When we tested the effects of 62 on the calculation of U2, 
we created a way in which one could compare a UI2 measurement 'corrected' for 62 to a model 
U2 and thus test whether or not a particular value of 62 was a good fit. In principle, one 
could calculate 0)2, varying the input value of 62, to find a best-fit value of 62- This best fit 
value of 62 would be independent of ag, as 62 affects the shape of 0)2, and ag affects only the 
normalization. Such an effort would require significant computational resources, and would 
certainly be prohibitive on the twenty-four separate subsamples we studied. Endeavoring to 
such a task using a sample that probes large volumes with minimal shot noise and precise 
photometric redshift estimations, such as the luminous red galaxies (LRGs), is more feasible 
(Ross et. al. 2007, in preparation). 
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7.2. C2 vs. bi 



Recent studies have found a relati onship between C 2 and b^ . iGaztanaga et al.l (120051 ) 
suggested the relationship C2 = &i — 1.2. iNishimichi et al.l (120061 ) measured results consistent 
with this relationship, and they showed that C2's dependence on bi is physically motivated. 
The major difference between these studies and ours is that we use galaxies over a much 
larger redshift range. Thus, we must determine if the theoretical relationship between C2 and 
bi is robust against changes in redshift. 



Following the methods of INishimichi et al.l (120061 ) . it is possible to calculate the first- 
and second-order bias parameters using a simple HOD model where the mean number of 
galaxies in a halo of mass M is given by 



(iV(M)) 



1 + 




M 
Ml 



M > M„ 
M < M„ 



(15) 



where a, Mi, and Mmin are free parameters and in general are fit such that this HOD model 
reproduces the measured clustering. Given this HOD model, the bias parameter bjsf is given 
by 



/^-^ dMn^UM, z) {N{M)) Bn{M, z) 



JN 



M (16) 
/^j;rfMn,„,,(M,z) (iV(M)) 

where nhaio is the mass fur iction of halos wit h mass m at redshift z determined via an ellip- 
soidal collapse model (e.g.. lSheth et al.ll200ll ) and Bn{m, z) is the nth-order bias coefficient of 



halos. We calculate both factors following the methods described in detail in INishimichi et al. 
(120061 ). but with cxg = 0.8 so that we can compare with the calculated bias with our bias 
measu rements (made assuming ag = 0.8). We adopt the fit parameters of IZehavi et al. 
J2005h determined for Mr < -19.5 (logioM^i^ = 11.76, logioMi = 13.15, a = 1.13) and 
Mr < -20.5 (logioM^i^ = 12.30, logioMi = 13.67, a = 1.21) for our HOD models which we 
use to calculate bi and C2 at the median redshifts of the different subsamples. 



In Figure [19] we display our calculated and measured C2 and bi , where the open green 
circles are the calculated values, black solid (all-), red dotted (early-), and blue dashed (late- 
type) crosses denote all of the measured values, and t he three lines repres ent C2 = 6i — 1.4, 
C2 = bi — 1.6, and C2 = 6i — 1.8. Using the HOD fits of IZehavi et al.l (120051 ) does not predict 
a strong relationship between C2 and bi. Changing the redshift weakly affects the calculated 
bias parameters, as the results show a significant change only when the HOD parameters are 
changed (the three calculations for Mr < —19.5 are hardly distinguishable, as are the three 
calculations for Mr < —20.5). Comparing these results to each of our measured C2 values for 
all galaxy types, the calculated values are consistent, though they are all greater than the 
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measured values. (The agreement is stronger if one considers that changing cxs will shift the 
values around.) Both the late-and all- type galaxies appear consistent with C2 being linearly 
dependent on bi, but the early- type galaxies C2 measurements show no clear dependence 
on bi. Our basic modeling suggests that there is nothing unusual about the relationship 
between C2 and 61 for early-type galaxies. It further suggests that for each galaxy-type there 
is redshift evolution of the HOD for a given luminosity. 



7.3. Late-type Galaxies 

The measurements made for late-type galaxies are significantly different for z < 0.3 
than for z > 0.3 (e.g., Figure [T3|) . The rise in correlation amplitudes at small scales, first 
reported by RossOG, happens only for z > 0.3. This suggests that late-type galaxies become 
much more likely to exist in close groupings at redshifts greater than 0.3. We are essentially 
measuring a preponderance of star forming galaxies in tight configurations as the redshift 
grows larger than 0.3. This in turn suggests that merger driven star formation becomes 
common at z > 0.3. This hypothesis is further supported by the fact that the bluest galaxies 
(L2) at z > 0.3 display the largest amplitudes at small scales (see Figure [H]), meaning that 
the galaxies with the most star formation are most likely to be found in tight groupings. 



This picture is broadly consistent with the concept of downsizing (e.g., ICowie et al. 



19961 ). which essentially states that higher-mass galaxies form stars earlier and more quickly 
than lower mass galaxies. The galaxies included in the volume limited samples are relatively 
bright {Mr < —19.5), and thus we find evidence of a preponderance of mer ger driven star 



formation only for z > 0.3. Our interpretation is supported by the results of iHeavens et al. 



(120041 ) who found that star formation peaked between a redshift of 0.3 and 0.8, using the 
fossil records of local SDSS galaxies. Further, they found the star formation rate to be a 
strong function of galaxy mass, implying that, for L > L*, there is little star formation at low 
redshift (this is supported by the dearth of late-type galaxies in the volume limited samples 
at z < 0.2 in our current analysis). This interpretation is further supported by recent results 



from the Cosmic Evolution Survey (jScoville et al. II2006I ). which found dramatically smaller 



star formation rates for galaxies 0.20 < z < 0.43 than for galaxies 0.43 < z < 0.65. 

The implications of C2jate being significantly smaller than that of early-type galaxies can 
be explained physically (see Tables 1 and 2). As the over-density of dark matter increases, 
the over-density of late-type galaxies becomes a smaller percentage of the dark mat ter over- 



densi ty. This is not a surprise, as the well-known morphology-density relationship (jPressler 



1980l ) tells us that the centers of clusters (i.e., the most over-dense regions) are filled with 



a smaller percentage of late-type galaxies than the outskirts of the clusters. It naturally 
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follows that late-type galaxies will have a smaller value of C2 than early-type galaxies. In 
terms of a halo occupation distribution (HOD), one would expect the fraction of late-type 
galaxies to decrease wi th the mass o f the h ost dark matter halo, which is the general trend 
recently determined by IZehavi et al.l (120051 ) . 



For z > 0.3, C2,iate is much closer to C2,eariy than any of the measurements made at 
lower redshift (see Tables 1 and 2). This suggests that at higher redshift, the fraction of 
red galaxies as a fu nction of d e nsity should display a shallower slope than at low redshift. 
This is observed by lYee et al.l (120051 ). as for galaxies with < —19.5, the slope in this 
relationship is on average smaller for galaxies with 0.4 < z < 0.6 than for galaxies with 
0.2 < z < 0.4. While these redshift ranges are different than the ones we employ, it confirms 
that the fraction of red galaxies versus density relationship shows a decre ase in slope as 



the re dshift increases (for < —19.5). This is in line with the results of iDressler et al. 
( 1l997l ). who found that at redshifts ~ 0.5 the fraction of spiral galaxies in clusters is 2-3 



times larger than in local clusters and that the spirals at higher redshift essentially replace 
the SO fraction. It is thus likely that the increase in C2,iate z > 0.3 is due to cluster spirals 
that have yet to evolve into SO galaxies. At smaller redshifts, the spirals have likely evolved 
into SO galaxies. Galaxies that we classify as late-type are therefore unlikely to be found in 
dense environments at low redshift, and thus C2^iate is significantly smaller for z < 0.3 than 
for z > 0.3. 



8. Conclusions 

The results presented in this paper represent the most complete and accurate determina- 
tion of the iV-th order correlations of photometrically selected galaxies. The measurements 
and the theoretical modeling used to interpret the measurements represent a significant im- 
provement over the Ross06 measurements. Taking all SDSS galaxies with 18 < r < 21 and 
measuring uj^ produces extremely interesting results, but it is only through volume limiting 
the sample and splitting by type, redshift, and luminosity that we are able to analyze the 
subtle effects that produce the measurements displayed for all galaxies. In doing so, we are 
able to quantify the nature of linear and non-linear clustering, and its dependence on type, 
redshift, and luminosity. 

We find that the linear bias parameter bi is smaller for late-type galaxies than for early- 
type galaxies, a result that is robust against changes in redshift and luminosity, but the ratio 
of bi^eariy to 6i,/ate docs Vary between 1.2 and 1.5 depending on the specific redshift /luminosity 
bin. We confirm that hi increases proportional to luminosity, as found in many previous 
studies. Significant evolution appears to occur in galaxies between a redshift of 0.3 and 0.4 
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as there is a large increase in bi going between galaxies with z < 0.3 and 0.3 < z < 0.4 but 
no significant change in bi between z = 0.2 and z = 0.3. 

The second-order bias, parameterized by C2, is significantly smaller for late- than early- 
types, and this is robust against any changes in the luminosity, redshift, or erg. This relation- 
ship between the non-linear bias of early- and late-type galaxies can be seen as a rigorous 
statistical restatement of the density - morp hology relationship and agrees with the results 



of the HOD analysis by IZehavi et al.l (120051 ). By applying a basic HOD model , we find our 



measured results are in fair agreement with the HOD parameters determined by lZehavi et al. 



fl2005[ ) 



There are large differences in the correlation measurements of late-type galaxies at 
redshifts greater- and less-than 0.3. This is broadly consistent with cosmic downsizing 



( ICowie et al.lll996l ). These differences predict a great amount of merger driven star formation 
at z > 0.3, and are consistent with the observed evolution in the density /morphology rela- 
tionship with redshift. Our results suggest that a detailed study of the density /morphology 
relationship as a function of redshift would find significant evolution at z ~ 0.3. 

If we require that bias be linear and set C2,au = 0, we find erg = 0.64]';q;q3, consistent with 
the lower limit on WMAP3 measurements constraine d by infiationary models. If instead, we 



set C2,aii equal to the cxg independent value found by iGaztanaga et al.l (120051 ) . we find that 
(Tg is a great match to the WMAP3-SDSS combined measurement. Considering all of the 
results, the most likely conclusion is that C2 is at least slightly negative and that erg < 0.8. 

As per usual, the results of our study demand more investigation. To this end, we are 
currently working to extend the analyses presented herein in two complementary directions. 
First, as was discussed earlier, the value of C2 can be constrained further by correcting the 
U2 measurements for the assumed C2 and measured bi/a^. Thus, by performing this analysis 
on the large, homogenous, photometric Luminous Red Galaxy sample from the SDSS and 
applying the results to our current analysis, we are thereby improving our measurements 
of higher-order bias terms and their dependence on galaxy- type, redshift, and luminosity. 
Second, we are also improving our theoretical interpretation of these results by performing 
a more rigorous halo occupation distribution model analysis of our higher-order correlation 
function measurements. 
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Table 1. The measured values of bias parameters bi and C2 for all measurements made 

using the fifteen z3 volume-limited samples. 



Type 


Mr Range 


z Range 


h ^Mpc Range 




bi 






C2 


XVDOF 


All 


< -19.5 


< 0.3 


8-36 


1.06 


± 


0.01 


1.10 


-0.45± 0.13 


0.44 


All 


-19.5 to -20.5 


< 0.3 


8-36 


0.97 


± 


0.01 


0.36 


-0.58± 0.20 


0.29 


All 


-20.5 to -21.5 


< 0.3 


8-36 


1.18 


± 


0.01 


0.76 


-0.38± 0.14 


0.16 


All 


< -19.5 


< 0.2 


8-24 


1.03 


± 


0.03 


0.24 


-0.36± 0.19 


0.52 


All 


< -19.5 


0.2 to 0.3 


8.5 - 38 


1.03 


± 


0.02 


1.13 


-0.32± 0.15 


0.29 


Early 


< -19.5 


< 0.3 


8-36 


1.33 


± 


0.03 


0.31 


0.06± 0.12 


2.40 


Early 


-19.5 to -20.5 


< 0.3 


8-36 


1.21 


± 


0.02 


0.51 


0.15± 0.15 


3.94 


Early 


-20.5 to -21.5 


< 0.3 


8-36 


1.50 


± 


0.03 


0.17 


0.08± 0.13 


0.34 


Early 


< -19.5 


< 0.2 


8-24 


1.27 


± 


0.04 


0.43 


0.10± 0.08 


1.23 


Early 


< -19.5 


0.2 to 0.3 


8.5 - 38 


1.27 


± 


0.03 


0.18 


0.05± 0.12 


1.21 


Late 


< -19.5 


< 0.3 


8-36 


0.87 


± 


0.01 


0.98 


-0.93± 0.34 


0.80 


Late 


-19.5 to -20.5 


< 0.3 


8-36 


0.84 


± 


0.01 


0.88 


-0.98± 0.27 


0.46 


Late 


-20.5 to -21.5 


< 0.3 


8-36 


1.04 


± 


0.02 


1.14 


-0.84± 0.25 


0.04 


Late 


< -19.5 


< 0.2 


8-24 


0.79 


± 


0.02 


0.10 


-1.03± 0.35 


0.27 


Late 


< -19.5 


0.2 to 0.3 


8.5 - 38 


0.84 


± 


0.01 


1.90 


-1.03± 0.35 


0.20 
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Table 2. The measured values of bias parameters bi and C2 for all measurements made 

using the nine zA volume-limited samples. 



Type 


Mj. Range 


z Range 


h ^Mpc Range 


bi 


XVDOF 


C2 


xVDOF 


All 


< 


-20.5 


< 0.4 


10 - 45 


1.35 ± 0.01 


1.73 


-0.29± 0.17 


0.89 


All 


< 


-20.5 


< 0.3 


8-36 


1.22 ± 0.03 


1.36 


-0.34± 0.21 


0.20 


All 


< 


-20.5 


0.3 to 0.4 


8-51 


1.39 ± 0.01 


2.81 


-0.17± 0.10 


0.46 


Early 


< 


-20.5 


< 0.4 


10-45 


1.64 ± 0.03 


0.31 


0.01± 0.14 


1.77 


Early 


< 


-20.5 


< 0.3 


8-36 


1.53 ± 0.03 


0.16 


0.08± 0.12 


0.70 


Early 


< 


-20.5 


0.3 to 0.4 


8-51 


1.66 ± 0.02 


1.12 


0.01± 0.08 


0.88 


Late 


< 


-20.5 


< 0.4 


10 - 45 


1.23 ± 0.02 


2.20 


-0.52± 0.27 


0.44 


Late 


< 


-20.5 


< 0.3 


8-36 


1.04 ± 0.02 


1.7 


-0.83d= 0.21 


0.03 


Late 


< 


-20.5 


0.3 to 0.4 


8-51 


1.25 ± 0.02 


3.66 


-0.38± 0.30 


0.37 
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0.2 0.4 0.6 



z 

Fig. 1. — Absolute r-band magnitude, M^, versus redshift, for galaxies in the DR5 PhotoZ 
catalog. The well-defined edge allows for the simple creation of volume limited samples. The 
magnitude-redshift areas of the two volume limited samples analyzed herein, z2> and zA, are 
outlined in red. 
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0.2 0.4 0.6 0.8 

z 

Fig. 2. — Normalized number of galaxies, n{z), in redshift bins Az — 10~^ determined by 
using the photometric redshifts of galaxies with 18 < r < 21 in the DR5 Photoz Table. This 
curve defines the redshift selection function we use to transform between real-space theory 
and our projected measurements. Similar curves are constructed for all other galaxy samples 
we analyzed herein. 
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Fig. 3. — The arca-averaged, angular, N-point correlation functions for 18 < r < 21 galaxies 

for N = 2, ,7 (left to right, top to bottom). Note the logarithmic scaling. Approximate 

power-law behavior is observed for all N to about 0.3 degrees, above which, however a power 
law clearly is no longer valid. 
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^(degrees) ^(degrees) 

Fig. 4. — The hierarchical amphtudes S3 through sj for galaxies with 18 < r < 19 (top 
left), 19 < r < 20 (bottom left), 20 < r < 21 (top right), and 18 < r < 21 (bottom right). 
The 6 values for each amplitude have been shifted and data with extremely large errors are 
omitted for clarity. The amplitudes are roughly constant, consistent with the hierarchical 
model, however significant structure is present. 
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Fig. 5. — The area-averaged angular N-point correlation functions for 18 < r < 21 early-type 



(red, t riangles) and late-type (blue, squares) galaxies (types determined by the lStrateva et al. 



(I2OOII ) u — r color cut), for = 2, ,7 (left to right, top to bottom). Note the logarithmic 

scaling. Clear differences between the two galaxy types are evident in both the slopes and 
amplitudes of these measurements. 
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Fig. 6. — The hierarchical amphtudes 5.s through 57 fo r late- type (left) and early-type (right) 
galaxies, with types determined by the lStrateva et al.l (I2OOII ) u — r color cut. The 6 values for 
S4 and Sj have been shifted and data with extremely large errors are omitted for clarity. At 
large scales, the early-types have significantly larger amplitudes than the late-types, implying 
a significant difference in bias between the two samples. 
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Fig. 7. — The hierarchical amphtude S3 for late-type (bottom) and early-type (top) galax- 
ies. T he solid line shows th e measurement made by using galaxies with types determined 
by the IStrateva et al.l (I2OOII ) color cut, while the points display galaxies whose types were 
determined using the type values from the DR5 Photoz table. 
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Fig. 8. — The hierarchical amphtudes S3 through S7 for all-type (left), late-type (middle) 
and early- type (right) galaxies for the volume limited sample of galaxies < z < 0.3, 
Mr < —19.5. The late-type galaxies show significantly different behavior than the late-type 
measurement made using the full sample of galaxies. The 6 values for S4 and Sy have been 
shifted slightly and data with extremely large errors are omitted for clarity. 
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Fig. 9. — The hierarchical amphtudes S3 (black triangles) and S4 (blue squares) for galaxies 
with < 2 < 0.3 split into -21.5 < Mr < -20.5 (left) and -20.5 < < -19.5 (right) 
displayed for all- (bottom), late- (middle), and early-type (top) galaxies. The 9 values for 
S4 have been shifted and data with extremely large errors are omitted for clarity (we also 
omit amplitudes S5 through Sj to be concise). The shapes are extremely similar, though 
the fainter galaxies display larger amplitudes, which is consistent with the more luminous 
galaxies having larger bias. 



-37- 



z<0.2 



0.2<z<0.3 



102 ^ I ■ 



10 



1 

102 



10 



1 

102 



10 



"T-m 1 1 — 1 I I I I 1 1 r 



I I I I I I II h 



H] — I I h I 



- I I II 1 1 1 1 1 — I I — ^ 



II 



H h 



4+1 



'"••'iiil 



*** 



II 



■■■ 



tl 



I ■ ■ 



102 



10 



102 



10 



102 



10 



0.1 1 

^(degrees) 



0.1 1 

^(degrees) 



10 



Fig. 10. — The hierarchical amphtudes S3 (black triangles) and S4 (blue squares) for galaxies 
with Mr < —19.5 separated into two redshift bins, < z < 0.2 (left) and 0.2 < z < 0.3 
(right) displayed for all- (bottom), late- (middle), and early- type (top) galaxies. The 9 values 
for S4 have been shifted and data with extremely large errors are omitted for clarity (we also 
omit amplitudes S5 through sy to be concise). 
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Fig. 11. — The average S3 measured from ten separate zS-like samples with redshifts de- 
termined by samphng the PDF of each galaxy's photometric redshift. The bold black error 
bars represent the standard deviation of the ten measurements, while the light red error-bars 
represent the jackknife errors of the S3 measurements using the true sample. 
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Fig. 12. — The hierarchical amphtudes S3 through sy for all-type (left), late-type (middle) 
and early-type (right) galaxies for the volume limited sample of galaxies < 2; < 0.4, and 
Mr < —21.5. These measurements look quite similar to those made on the full sample. The 
6 values for S4 and S7 have been shifted slightly and data with extremely large errors are 
omitted for clarity. 
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Fig. 13. — The hierarchical amphtudes S3 (black triangles) and S4 (red squares) for galaxies 
with Mj. < —20.5 separated into two redshift bins, < z < 0.3 (left) and 0.3 < 2; < 0.4 
(right) displayed for all- (bottom), late- (middle), and early- type (top) galaxies. The 9 values 
for S4 have been shifted slightly and data with extremely large errors are omitted for clarity 
(we also omit amplitudes S5 through S7 to be concise). 
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Fig. 14. — The hierarchical amphtudes S3 (black triangles) and S4 (blue squares) for galaxies 
with Mr < —20.5 and 0.3 < z < 0.4 displayed for Ll (closed symbols) and L2 (open symbols) 
galaxies. The 6 values for S4 have been shifted slightly and data with extremely large errors 
are omitted for clarity. 
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Fig. 15. — The S3 measurements for Early-type galaxies from the z < 0.2, Mr < —19.5 
(black triangles) and the 0.3 < z < 0.4, Mr < —20.5 (red circles) samples plotted against 
the angular scale (left panel) and the equivalent physical scale (right panel), with the lower 
redshift galaxies' S3 divided by the relative bias of the two samples (1.31). When plotted 
on the physical scale, the shape of the measurements agree, indicating that the features 
measured are physical in nature, and not a systematic effect. 
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Fig. 16. — Same as Figure [13] only early- and late-type galaxies are split by the color method 
instead of the photo method. This suggests that the color method does not adequately 
separate early- and late-type galaxies for z > 0.3. 
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Fig. 17. — The ratio of uj2,c to uj2-, where c<)2,c calculated by correcting over-densities for a 
second-order bias term (62) equal to -0.3 and UJ2 is calculated in the standard way. Both 
are calculated using the zi sample of galaxies with —20.5 < < —19.5. The ratio begins 
to grow significantly greater than 1 for ^ < 0.66°, corresponding to a physical scale of 
approximately 8/i^^Mpc. This corresponds with the the lower bound of the weakly non- 
hnear regime. All bias values are thus calculated for r > 8/i~"^Mpc. 
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Fig. 18. — The 1 cr allowed regions of parameter space for cxg and for all- (black, solid), 
early- (red, dashed), and late-type (blue, dotted) galaxies. The allowed second-order bias 
varies greatly with erg, but the differences in C2 as a function of galaxy type are robust. 
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Fig. 19. — The measured for each subsample pfotted against the measured h\. Each point 
is marked by it error-bars with early-type galaxies being red and dotted, late-type galaxies 
blue and dashed, and all-type galaxies being black and solid. The solid black lines represent 
C2 = &i — 1-4, C2 = &i — 1.6 and C2 = h\ — 1.8. The g reen open cir c les di splay the calculated 
C2 and h\ using the HOD parameters determined by IZehavi et al.l (120051 ). 



